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Investigation of a-induced reactions on 127 I for the astrophysical 7 process 
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BACKGROUND: The 7 process in core-collapse supernova explosions is thought to explain the 
origin of proton-rich isotopes between Se and Hg, the so-called p nuclei. The majority of the reaction 
rates for 7-process reaction network studies has to be predicted in Hauser-Feshbach statistical model 
calculations. Recent investigations showed problems in the prediction of a widths at astrophysical 
energies. This impacts the reliability of abundance predictions in the upper mass range of the p 
nuclei. 

PURPOSE: Measurement of the 127 I(q,7) 131 I and 127 I(a,n) 130 I reaction cross sections close to 
the astrophysically relevant energy range to test the predictions, to derive an improved reaction 
rate, and to extend the database required to define an improved global optical a+nucleus potential. 

METHODS: The cross sections were derived using the activation technique, the yield of the 
emitted 7 and characteristic X-ray photons were measured using a LEPS and a HPGe detector. 

RESULTS: The cross sections of the 127 I(q,7) 131 Cs reaction have been determined for the first 
time, at energies 9.50 < E c . m . < 15.15 MeV. The 127 I(a,n) 130 Cs reaction was studied in the range 
9.62 < E c .m. < 15.15 MeV. Furthermore, the relative intensity of the 536.1 keV 7 transition was 
measured precisely, its uncertainty was reduced from 13% to 4%. The results were then compared 
to Hauser-Feshbach calculations which were also used to extend the cross sections into the astro- 
physically relevant region and to compute the reaction rate. 

CONCLUSIONS: The comparison to statistical Hauser-Feshbach model calculations showed 
that the a width can be described well in the measured energy range using a standard, energy- 
independent global optical potential. The newly derived stellar reaction rates at 7-process tempera- 
tures for 127 I(a,7) 131 I and its reverse reactions, nevertheless, are faster by factors 4 — 10 than those 
from previous calculations, due to further improvements in the reaction model. The importance 
of the inclusion of complete level schemes into the Hauser-Feshbach calculations is illustrated by 
comparing the impact of two different level schemes, one of them extending to higher excitation 
energy but not containing all relevant levels. 

PACS numbers: 26.3Q.-k, 25.55.-e, 27.60.+j 
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I. ASTROPHYSICAL MOTIVATION 

The bulk of the stable isotopes heavier than iron ob- 
served in our Solar System are produced via neutron cap- 
ture reactions (IJ. In order to reproduce the abundance 
pattern of the heavy isotopes, two different neutron cap- 
ture processes have to be assumed: the s-process (which 
has a main and a weak component [3]) and the re- 
process @, m. 

On the proton-rich side of the valley of stability, how- 
ever, there are 35 nuclei which are separated from the 
path of the neutron capture processes. These mostly 
even-even isotopes between 74 Se and 196 Hg are the so- 
called p nuclei [6| . It is generally accepted that the main 
stellar mechanism synthesizing the p nuclei - the so- 
called 7 process - involves mainly photodisintegrations, 
dominantly (7,11) reactions on preexisting, neutron-richer 
s and r seed nuclei. The high energy photons - necessary 
for the 7-induced reactions - are available in explosive 
nucleosynthctic scenarios like a core-collapse supernova 



* ggkiss@atomki.hu 



Shockwave moving through the Ne/O layer of a massive 
star and reaching temperatures around a few GK. Favor- 
able conditions have recently been found also in Type la 
Supernovae Q. 

Consecutive (7,11) reactions can drive the material to- 
wards the proton rich side of the valley of stability. As 
the neutron separation energy increases along this path, 
(7,p) and (7,0) reactions become stronger and bend the 
reaction flow towards lighter elements [a-[lO{. Theoreti- 
cal investigations agree that (7,p) reactions are more im- 
portant for the lighter p nuclei while (7, a) reactions are 
mainly important at higher masses (A > 100). However, 
calculations based solely on the 7 process in massive stars 
are unable to reproduce the observed abundance pattern 
of the p nuclei. To solve this problem, several processes 
have been suggested to fill in the missing nuclear abun- 
dances, e.g., the rapid proton capture (rp) process on 
the surface of accreting neutron stars [II], [l2| and the 
neutrino-induced p process (yp process) in the deepest 
layers of a core-collapse supernova ejected in a neutrino- 
wind [ljj H3 |. Considerations regarding isotope ratios 
in meteoritic material, however, severely constrain their 
possible contribution 15 1. 
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Theoretical studies of the nuclear uncertainties in the 
nucleosynthesis of the p nuclei make use of large reaction 
networks with mainly theoretical reaction rates (taken 
from the Hauser-Feshbach (H-F) model [l6j]). They have 
shown that the reaction flow for the production of heavy 
p nuclei (140 < A < 200) is strongly sensitive to the 
(7, a) photodisintegration rates [H, [To] . Recent experi- 
ments, however, indicate that the H-F predictions may 
overestimate the a-capture cross sections at low energies 
by factors 3 — 10. This would also strongly impact the 
astrophysical reaction rates. Unfortunately, the avail- 
able data are scarce, only a handful of reactions have 
been studied at sufficiently low energies [l7l - [26j . Conse- 
quently, experimental data is urgently needed to confirm 
the path of the 7 process, especially for reactions involv- 
ing a particles. 

It has to be understood that straightforward measure- 
ments of (7, a) reactions on target nuclei in the ground 
state (g.s.) cannot provide the required information for 
astrophysical (7,0) reactions in a stellar plasma, involv- 
ing a large number of reactions on nuclei in thermally 
excited states. The sensitivity of the photodisintegra- 
tions on g.s. nuclei is very different from the one of stel- 
lar rates and measurements, therefore test the prediction 
of different nuclear properties than actually entering the 
calculation of the rates [13, HH . In a high-temperature 
stellar plasma, as encountered in the 7 process, there is 
complete detailed balance between forward and reverse 
rates, therefore (7, a) rates can be easily converted to 
(a, 7) ones. This is only possible with stellar rates, how- 
ever [27],[29|,[3(|. If the g.s. contribution to the stellar rate 
< X < 1 is close to unity, i.e., reactions on the g.s. of 
the target nucleus provide most of the rate, a measure- 
ment will be able to directly provide the rate [2^, l3lj |. 
Photodisintegrations have several orders of magnitude 
lower X than captures [28J, |32|. Therefore, a measure- 
ment of (a, 7) is required to test the predictions of stellar 
(7, a) rates. 

It is essential, however, to measure as close as possible 
to the astrophysically relevant energies, i.e., the energies 
from which the largest contribution to the reaction rate 
integral are coming [27|, |33[ . Because of the strong en- 
ergy dependence of charged-particle reaction widths, the 
dependences of the reaction cross sections at astrophys- 
ical, subCoulomb energy and those at higher, measured 
energy can be quite different. For example, all astrophys- 
ical a captures are determined by the a width, whereas 
well above the Coulomb barrier the 7 width is mostly 
dominating and other particle widths may play a role 
close to channel openings [13, HM|- If a measurement 
is not possible in the a dominated energy region, the 
simultaneous measurement of a reaction channel domi- 
nated by the a width may prove useful, e.g., the (a,n) 
channel. The problem with this is that the (a,n) reac- 
tion on proton-rich target nuclei, interesting for the 7 
process, always exhibits a negative Q value, thus limit- 
ing the measurable energy range again to higher energies 
than astrophysically relevant. Nevertheless, this can be 



used as an initial test of the predicted a width. As further 
discussed in Sec. IIII1 it remains important to remember 
that a discrepancy between theory and experiment in a 
capture does not necessarily imply an incorrect predic- 
tion of (a, 7) and (7,0?) reaction rates in stellar plasmas. 

In order to extend the experimental database for the 
astrophysical 7 process and to test the reliability of H-F 
statistical model predictions in the A > 100 mass range, 
the cross sections of the 127 I(a,n) and 127 I(a,7) reac- 
tions have been measured for the first time, using the 
activation method. The 127 I(a,n) reaction was studied 
in the range 9.62 < E c . m . < 15.15 MeV. The radia- 
tive a-capture cross sections were measured at 9.50 < 
E c . m . < 15.15 MeV. This is close to the astrophysically 
relevant energy region (the Gamow window), which cov- 
ers 5.8 — 8.2 MeV at a plasma temperature of T = 2.5 
GK g3. 

The paper is organized as follows: The experimental 
details are described in Sec. [TTJ the resulting cross sec- 
tions are presented in Sec. IIII1 where also a comparison 
to predictions is shown and the implications for the as- 
trophysical rate are discussed; Sec. IIVI provides the con- 
clusions and a short summary. 



II. EXPERIMENTAL TECHNIQUE 

The element iodine has only one stable isotope 127 I. 
The (01,11) and (a, 7) reactions studied in the present work 
lead to unstable cesium ( 130,131 Cs) isotopes. The 130 Cs 
reaction product has a half-life of 29.21 minutes and its 
electron capture decay is followed by the emission of a 
536.1 keV 7 ray. The half-life of the 131 Cs produced 
via the 127 I(a,7) reaction is 9.69 days. It decays exclu- 
sively by electron capture and no 7 radiation is emitted 
during its decay. This electron capture is followed by 
the emission of characteristic X-rays. For the determina- 
tion of the alpha-capture cross sections the yield of these 
characteristic X-rays was used similarly as in [lEJ. The 
technical details of this approach are described in [36| . 
Further experimental details are described below in de- 
tails. The decay properties of the reaction products are 
summarized in Table HI 



A. Target production and irradiation 

Because of its high electronegativity, iodine reacts vi- 
olently with alkali metals such as Na and K. Therefore, 
it naturally occurs mostly in compounds. For the target 
production a KI compound was used because it is chemi- 
cally stable and its melting point is sufficiently high. The 
targets were made by vacuum evaporation from a Mo 
boat onto 2 /im thick aluminum foils. During the evap- 
oration, the distance between the Al backing and the 
evaporation boat was 9 cm. Using this relatively large 
distance, uniform targets could be produced. 



3 



TABLE I. Decay parameters of the 130 ' 131 Cs reaction products, formed via the 127 I(a,n) 130 Cs and the 127 I(a,7) 131 Cs reactions. 
The relatively weak characteristic X-ray transitions (indicated by italics) were only used to measure the relative intensity 
of the 536.1 keV 7 transition. 



Product 
nucleus 


Decay 
mode 


Half-life 
(hour) 


transition 


X- and 7-ray 
energy (keV) 


Relative 
7-intensity 
per decay (%) 


Ref. 


131 Cs 


e 100% 


232.54 ± 0.38 




29.461 


21.1 ± 0.5 


[34] 








K q i 


29.782 


38.9 ± 0.9 




130 Cs 


e 96.4% 


0.4868 ± 0.00067 


k q2 


29.461 


21.6 ± 0.7 


f35] 








K q i 


29.782 


11.7 ± 0.6 












33.56 


2.05 ± 0.05 












33.62 


3.95 ± 0.10 












34.42 


1.20 ± 0.03 












536.1 


3.80 ± 0.50 














3.81 ± 0.15 


present work 



The absolute number of target atoms was determined 
by PIXE (Proton Induced X-ray Emission) method [37} 
at the PIXE chamber installed on the left 45° beamline 
of the 5 MV VdG accelerator of ATOMKI, where the 
current of the beam and the dead time can be measured 
very precisely. A more detailed description of the setup 
can be found in [38| . A homogeneous beam of 2 MeV 
protons with 0.5 mm diameter and « 1 nA current was 
used for the thickness measurement. The total collected 
charge in the case of each target was about 1 /iC. A typ- 
ical PIXE spectrum is shown in Figure 1. The spectrum 
was fitted using the PIXEKLM program code [39]. The 
peaks used for the analysis are marked. 

The thickness of the KI targets were between 106 
and 507 /^g/cm 2 , corresponding to about 5.3 x 10 17 
I atoms/cm 2 and 24.0 x 10 17 I atoms/cm 2 . The X- 
ray attenuation factor for the K a i.2 characteristic X-rays 
emitted during the decay of the alpha-capture reaction 
products was calculated using the LISE code and 
was found to be below 0.4% for these target thicknesses. 
Since the chemical composition of the target is impor- 
tant for the energy-loss calculations, the ratio of the K 
and I atoms in the target were also determined and were 
found to be 1:1 (within 1%). The precision of the deter- 
mination of the number of target atoms was better than 
4%. Using the PIXE method, the following impurities 
have been found in the target and the backing: Fe, Cu, 
Ga, Zn, V (less than 50 ppm each). To check the PIXE 
results, weighing was also used to determine the number 
of target atoms in a few cases. The agreement between 
the results of the two methods is within 5%. 

The KI targets have been irradiated with a beams from 
the MGC20 cyclotron of ATOMKI. The energies were in 
the range of 9.81 < E a < 15.5 MeV, covered in steps 
of about 0.5 MeV. After the beam-defining aperture, the 
chamber was insulated and a secondary electron suppres- 
sion voltage of —300 V was applied at the entrance of the 
chamber [221 ]. The typical beam current was between 0.1 
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FIG. 1. (Color online) PIXE spectrum measured by bom- 
barding the KI targets with 2 MeV protons. The peaks used 
for the analysis are marked. Peaks belonging to impurities in 
the target and/or the backing are indicated, too. 
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FIG. 2. Off-line 7 - (upper panel) and characteristic X-ray 
(lower panel) spectra normalized to the length of the count- 
ings, taken after irradiating a 127 I target with E a = 13.0 
MeV(E c . m . = 12.58 MeV) beam. The 7-line and X-ray tran- 
sitions used to measure the cross section of the 127 I(«,n) 130 Cs 
and 127 I(q, 7) 131 Cs reaction cross sections as well as the main 
background lines are marked. 



and 0.6 p/iA, the length of each irradiation was between 
0.92 and 24 hours, corresponding to about 2.7 x 10 15 and 
2.7 x 10 17 total incident a particles. Several beam tests 
were performed to check the target stability before the 
experiment. These tests showed that there was no dete- 
rioration of the targets using a beam current less than 
0.8 pfiA. The current integrator counts were recorded in 
multichannel scaling mode, stepping the channel every 
minute to take into account the changes in the beam 
current. 



B. 7- and characteristic X-ray countings 

For the 7 counting a 100% relative efficiency HPGe 
detector placed in ultra low background (ULB) configu- 
ration was used. The X-ray countings were carried out 
using a Low Energy Photon Spectrometer (LEPS). In 



order to reduce the laboratory background a multilayer 
quasi 47T shield has been built around this detector in- 
cluding an inner 4 mm thick layer of copper and a 2 mm 
thick layer of cadmium and an 8 cm thick outer lead 
shield |41| . The low yields measured in the present work 
necessitated the use of small source-to-detector distances 
both for the 7- and X-ray countings. The distance be- 
tween the target and the endcup of the ULB detector was 
1 cm. In the case of the X-ray countings, the targets were 
placed 3 cm from the Be entrance window of the LEPS 
detector. 

The detector efficiencies for both detectors used in the 
present work had to be known in close geometry with 
high precision. A similar procedure had been used to de- 
rive the photopeak efficiencies in the case of the HPGe 
and LEPS detectors. At first the absolute detector effi- 
ciency was measured in far geometry: at 10 cm distance 
from the surface of the detector, using calibrated 57 Co, 
60 Co, 133 Ba, 137 Cs, 152 Eu, and 241 Am sources with the 
HPGe detector, and 133 Ba, 152 Eu, and 241 Am sources 
with the LEPS detector. Since the calibration sources 
(especially 133 Ba, 152 Eu) emit multiple 7-radiations from 
cascade transitions, in the close geometry a strong true 
coincidence summing effect is expected resulting in an 
increased uncertainty of the measured efficiency. There- 
fore, no direct efficiency measurement in close geometry 
has been carried out. Instead, the activity of several ir- 
radiated 127 I targets has been measured both in close 
and far geometry. Taking into account the time elapsed 
between the two countings, a conversion factor of the effi- 
ciencies between the two geometries could be determined 
and used henceforward in the analysis. 



1. Measurement of the 127 1(a,n) 130 Cs reaction products 

The relative intensity of the 536.1 keV 7 transition - 
emitted during the decay of the 130 Cs reaction product - 
is known with limited precision, its uncertainty is 13%. 
In order to increase the precision of the cross section data 
derived in the present work, first, the uncertainty of the 
branching ratio of this transition was reduced. Two rel- 
atively thin KI samples have been irradiated with a E a 
= 15 MeV (E c , m . = 14.52 MeV) beam and their spectra 
were measured using the LEPS detector. The yield of 
the 536.1 keV 7 transition was normalized to the yield 
of the K a i j2 , to the yield of the K Q i j3 and to the yield of 
the K^2 X-ray transitions, taking into account their well 
known relative intensities and the detector efficiencies. 
The weighed average of the 6 normalized relative inten- 
sities was found to be 3.81% ± 0.15%. Its uncertainty 
was calculated by accounting for the following partial un- 
certainties: uncertainty of the detector efficiency ratios 
(3%), uncertainty of the X-ray relative intensities taken 
from the literature (< 0.5%) and statistical uncertainty 
(< 2.0%). 

The activity of the 127 I(a,n) reaction products was 
measured after the end of the irradiation by counting the 
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TABLE II. Measured cross sections of the 127 I(a,n) 130 Cs re- 
action. 



E 


cross section 


(MeV) 


(mbarn) 


9 f)2 i 10 


013 + 002 


mm +nn<) 


045 ± 007 


10 64 ± 06 


17 ±0 01 

V / ■ 1 I — 1 — V 7 a\J I 


1 71 + 08 

-L\J . 1 _L 1 VJ.yjlJ} 


1 8 ± 05 


11.10 ±0.09 


0.59 ±0.05 


11 13 ± 06 

1 1.1'/ _i_ UiUU 


66 ± 11 


11.58 ±0.06 


1.61 ±0.12 


11 60 ± 06 


1 53 ± 27 

1 ■ ' I * } — 1 — \J a £u 1 


12.18 ± 0.09 


4 30 ± 3 


1 9 99 + n 06 


5 01 + 89 


12.58 ±0.07 


8.78 ±0.61 


12.59 ±0.09 


10.2 ±1.6 


13.18 ±0.09 


17.1 ±2.5 


13.55 ±0.07 


32.1 ±4.8 


14.19 ±0.09 


68.6 ±10.1 


14.52 ±0.07 


119.9 ±8.4 


14.53 ±0.08 


115.6 ±17.0 


15.15 ±0.10 


167.6 ±24.7 



yield of the 536.1 keV 7 transition at and above E c . m . 
= 10.71 MeV. The upper part of Fig. [5] shows a typical 
spectrum collected using the HPGe detector after the ir- 
radiation of a KI target with 13 MeV a beam (-E c .m. = 
12.58 MeV). At the low energy irradiations (below E c m . 
= 10.71), the yield of the 536.1 keV transition was not 
sufficient for the analysis. For this reason, the emitted 
Kai,2 characteristic X-ray photons were counted. To de- 
termine the cross sections at these energies, the length 
of the irradiations was short (about 1 hour), in order 
to limit the beam-induced background. To verify this 
approach, several irradiations at higher energies were re- 
peated, the 127 I(a,n) reaction cross sections were mea- 
sured two times, first by counting the yield of the 536.1 
keV transition, and in another irradiation by measuring 
the yield of the emitted characteristic K Q i :2 X-ray pho- 
tons. The experimental cross sections are given in Table 
HIl Between the end of the irradiation and the 7 / X-ray 
counting, a waiting time of about 0.25 hours was inserted 
in order to decrease the yield of the disturbing short-lived 
activities. All spectra were taken for two hours at most 
and stored regularly in order to follow the decay of the 
130 q s reac ^j on product. 



2. Measurement of the 127 1(a^) 1Z1 Cs reaction products 

The yield of the characteristic X-rays — emitted dur- 
ing the electron-capture decay of the produced unstable 
cesium isotopes — were used for the determination of the 
(a, 7) and (ev,n) reaction cross sections. The disadvan- 
tage of this approach is that it is not able to distinguish 
between the decay of the different isotopes of the same 
element. However, using a similar approach as in [25| . 



TABLE III. Measured cross sections of the 127 I(a, 7) 131 Cs 
reaction. 



E 


cross section 


(MeV) 


(mbarn) 


9 50 ± 12 


18 ± 04 


10 00 ± 09 


65 ± 12 


10 71 ± 08 


1 76 ± 29 


11 1 q ± n dfi 

_i_j_._i_o zn u.uu 


qn _l n 7Q 
o.yu m U- / y 


11.60 ±0.06 


8.54 ±0.74 


12.22 ±0.06 


21.6 ±1.7 


12.59 ±0.09 


33.4 ±2.6 


13.18 ±0.09 


42.8 ±3.5 


13.55 ±0.07 


63.5 ±5.0 


14.19 ±0.09 


96.6 ±7.3 


14.53 ±0.08 


123.2 ±9.2 


15.15 ±0.10 


127.7± 11.0 



the two open reaction channels can be separated. Cal- 
culations show that close above the threshold the (a,n) 
channel becomes dominant, its cross section exceeds by 
orders of magnitude that of the (a, 7) channel. Since the 
half- life of 131 Cs is about 478 times longer than that of 
130 Cs, if the X-ray countings are carried out more than 
1 day after the end of the irradiation, the target's 130 Cs 
activity decreases to a negligible level and the decay of 
the longer- lived 131 Cs can be measured. 

The energies of the emitted characteristic X-ray K Q i.2 
lines are 29.782 and 29.461 keV, respectively. Since the 
resolution of the LEPS detector is typically between 400 
eV (for a 5.9 keV 7 line) and 680 eV (for a 122 keV 7 
ray), in the X-ray spectra it is not possible to distinguish 
between K a i and K a 2 transitions. Instead, the sum of the 
emitted characteristic X-rays was used for the analysis. 
The X-ray spectra were taken for 1-6 days and stored 
regularly in order to follow the decay of the 131 Cs reaction 
products. In most of the cases, the X-ray counting was 
carried out twice and consistent results were found. The 
lower part of Fig. [2] shows a typical spectrum collected 
using the LEPS detector after the irradiation of a KI 
target with a E a = 13 MeV {E c . m . = 12.58 MeV) beam. 

The experimental cross sections for this reaction chan- 
nel are given in Table IIIII 



III. RESULTS AND DISCUSSION 

The measured ev-induced cross section values are listed 
in Tables ITT1 and Hill The quoted uncertainty in the E C An , 
values corresponds to the energy stability of the a-beam 
and to the uncertainty of the energy loss in the target, 
which was calculated using the SRIM code [13] • Some ir- 
radiations were repeated at the same energies. The cross 
sections are then derived from the averaged results of 
the irradiations weighted by the statistical uncertainty 
of the measured values. The uncertainty of the cross sec- 
tions is the quadratic sum of the following partial errors: 
efficiency of the HPGe detector and LEPS (6 and 5 %, 
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respectively), number of target atoms (5%), current mea- 
surement (3%), uncertainty of decay parameters (< 13.2 
%) and counting statistics (1.6 - 17.9%). 



A. Astrophysical implications 

In astrophysical investigations it is common to quote 
the astrophysical S factor S(E) which removes the energy 
dependence due to the projectile penetration through the 
Coulomb barrier from the cross sections <r(E) [43j |. 



27T7/ 



S(E) = Ea(E)e 
with r\ being the Sommerfeld parameter 



V 



= Z p Z T e 2 f fi \ 



1/2 



\2EJ 



(1) 



(2) 



The charge numbers Z p , Zt of projectile and target, re- 
spectively, and their reduced mass fx enter the Sommer- 
feld parameter. When there is an uncertainty in the en- 
ergy E, the calculation of the S factor errors becomes 
more complicated because the energy enters also via the 
Sommerfeld parameter. The error bars on a energy and 
cross section translate into an error region for the S fac- 
tor which is of trapezoid shape, instead of the usual rect- 
angular shape. These error trapezoids are also used in 
Figs. [3] and 31 showing a comparison of the experimental 
S factors for the (a,j) and (a,n) reactions, respectively, 
with theoretical calculations in the H-F approach. 

Figure [3] compares the standard NON-SMOKER pre- 
diction [44[ (dashed line) to the (a/y) data. The calcu- 
lation overestimates the experimental S factor by a fac- 
tor of two at the upper end of the investigated energy 
range and a factor of about 6 at the lower end. Also 
shown are calculations with the code SMARAGD [45j . 
version 0.8.4s, using different numbers of discrete excited 
states (this issue is discussed below in Sec. IIIIBj) . The 
SMARAGD calculation with 5 discrete states included is 
similar to the NON-SMOKER values, although slightly 
closer to the data. As pointed out in Sec. HI however, 
one should not jump to conclusions regarding the pre- 
diction of the stellar a capture rate. The sensitivity of 
the stellar reaction rate to a variation of the averaged 
7, neutron, and a widths by a factor of two is shown in 
Fig. [5j (All discussed cross sections and rates are insen- 
sitive to a variation of the proton width, therefore it is 
never shown.) A sensitivity value \s\ = 1 implies that the 
change in the cross section or rate is of the same mag- 
nitude as the variation factor, a value s — shows that 
there is no dependence on this width 28j. The sign of 
the plotted sensitivities informs whether the cross sec- 
tion or rate changes in the same direction as or opposite 
to the width. We note that the (a, 7) cross sections and 
rates change inversely proportional to the neutron width 
whereas they depend proportionally on the 7 and/or a 
width. The dependences are reversed between neutron 
and 7 width in the (a,n) channel. For a detailed expla- 
nation of the possible cases, see [28j |. 



As can be seen in Fig. O the reaction rate is only sen- 
sitive to the a width at the temperature region relevant 
in the 7 process. Comparing this to the sensitivity of the 
capture cross sections in Fig. [S] it is immediately obvi- 
ous that, in addition to the a widths, those also strongly 
depend on the 7 and neutron widths in the measured en- 
ergy range. Only below the (a,n) threshold, the situation 
is similar to the one of the rate. Therefore it is impos- 
sible to further disentangle the different contributions of 
the widths to the total deviation from experiment and to 
draw strong conclusions on the basis of the experimental 
capture data alone. 

The (a,n) cross sections (and S factors) at the upper 
end of the measured energy range, on the other hand, also 
are only sensitive to the a width, similar to the reaction 
rate, as seen in Fig. [71 Inspection of Fig. [4j comparing 
the (a,n) data to the calculations, finds that the theory 
values are close to the experimental ones, especially at 
higher energies. While NON-SMOKER is slightly above 
the data, the SMARAGD calculation using five excited 
states is in good agreement with the measured S factors. 
It is to be concluded, therefore, that the a width is pre- 
dicted well by the SMARAGD code and any deviations 
from the (a/y) experiment are due to problems in the 
neutron and/or 7 widths. Without further information, 
it is impossible to further identify one of these widths as 
main source of error. 

Since the g.s. contribution to the stellar rate is X = 
0.35 at 2.5 GK (28|, the experimental data cannot be 
used to compute the stellar capture and photodisinte- 
gration rates. The a width, however, was shown to be 
the important width for the rate and its prediction was 
confirmed by the data, at least in the measured energy 
range. Therefore the stellar rate can be calculated uti- 
lizing the same approach that was used for comparison 
with the data but including the thermal plasma effects. 
The SMARAGD code uses a superior method to calcu- 
late transmission coefficients and averaged widths below 
the Coulomb barrier compared to the NON-SMOKER 
code [274 and we used this code (with five discrete ex- 
cited states included) to provide the stellar reactivities in 
Table [TV] The REACLIB parameters (H for a fit based 
on the new reactivity values are given in Table |V] The 
fit accuracy is better than 0.2% within the relevant tem- 
perature range. 

Figure [5] shows a comparison of the previous standard 
rate from [44j (which was fitted to provide the REACLIB 
parameters in [46]). At 2.0 and 2.5 GK, the rate is in- 
creased by factors 10 and 4, respectively. This is neither 
a consequence of different optical a+nucleus potentials 
(since both codes use the same 47]) nor of a modifica- 
tion based on the new experimental data. It simply shows 
the difference in the numerical solution of the Schrodinger 
equation at strongly subCoulomb energies [271 ] . 

Finally, it has to be cautioned that the a width is still 
not well constrained at astrophysically relevant energies 
by the experiment. It is conceivable, for instance, that 
the subCoulomb energy dependence is different than the 
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FIG. 3. (Color online) Experimental S factors of the 127 I(a,7) 131 Cs reaction compared to theoretical predictions with the 
NON-SMOKER and SMARAGD codes (see text). 



one found in the measured energy range. In fact, there 
could be an indication of this at the lowest two data 
points. Inspecting Fig. 21 there seems to be a slight 
deviation below 11 MeV. One might be tempted to as- 
sign this to the same problem with the neutron and/or 
7 width showing up in the capture data because it ap- 
pears just as the sensitivity to the neutron and 7 widths 
increases. Comparing the sensitivities of the (a, 7) and 
(a,n) reaction, however, reveals that the required modifi- 
cation would be different in the two cases. The calculated 
(a, 7) excitation function could be shifted to lower values 
by either increasing the 7 width or decreasing the neu- 
tron width. Both changes would act oppositely in the 
(a,n) excitation function, moving the theory values fur- 
ther away from the data at low energy. Therefore, if this 
is the sign of a real effect, it has to be due to a changed 
energy dependence in the a width. The indication for 
this in the current data, however, is too vague to draw 
any conclusion. 



B. The problem of the selection of discrete excited 
states in Hauser-Feshbach calculations 

These reactions are also well suited to discuss another 
problem appearing in the prediction of reactions, es- 
pecially for unstable nuclei close to stability with par- 
tially known level schemes. The H-F transmission co- 
efficients, required to compute the total width in each 
reaction channel, are calculated from sums of individual 
transitions from a compound nuclear state to all ener- 
getically possible final states. In each reaction channel 
excited states have to be considered up to an energy 
E m ax = E c . m . + Q c . Depending on the reaction Q value 
of the channel Q c , often discrete excited levels are only 
known up to an energy i?i as t < £ max . Then the calcu- 
lation of transitions to levels above £"i as t are performed 
by an integration over a nuclear level density instead of 
a sum over discrete, isolated states [27l[28|. 

Far off stability, only few or no discrete states are 
known and the level density integration is used from low 
excitation energies onwards or directly above the g.s. But 
also close to or at stability, where levels are often speci- 
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FIG. 4. (Color online) Experimental S factors of the 127 I(a,n) 130 Cs reaction compared to theoretical predictions with the 
NON-SMOKER and SMARAGD codes (see text). 
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5. Sensitivity of the 127 I(a,7) 131 Cs stellar reaction rate 



variation of the 7, neutron, and a width as function 
of plasma temperature T. The sensitivity to a variation of 
the a width is constantly unity at all temperatures. The 
temperature relevant for this reaction in the astrophysical 7 
process is 2 < T < 2.5 GK. 



fied up to .Emax, it is necessary to choose an appropriate 
cut-off at an energy below E max because it is essential for 
a successful prediction to have a complete level scheme 
included. Often, even when levels are given up to high 
excitation energy in a database like NuDAT or ENSDF 
[H,[49j], these are only partial level schemes, derived from 
a limited number of experiments or coming from informa- 
tion on rotational bands. At higher excitation energies it 
becomes also harder to experimentally resolve individual 
levels lying close to each other. An educated guess has to 
be made when preparing the input for H-F calculations, 
estimating up to which excitation energy the experimen- 
tal level scheme available in databases can be consid- 
ered complete. This is especially important in particle 
channels because particle transitions to low-lying excited 
states contribute most, whereas in the 7 channel tran- 
sitions to higher excitation energies (with lower relative 
7 energy) are dominant [28l . [50| and lead to the region 
where the nuclear level density is used, anyway. 

The Q value of the reaction 127 I(a,n) 130 Cs is -7.729 
MeV and therefore the impact of the low-lying excited 
states is even more pronounced. Table IVTl lists the set of 
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FIG. 6. Sensitivity of the 127 I(a,7) 131 Cs reaction cross sec- 
tions to a variation of the 7, neutron, and a width as function 
of cm. energy. The sensitivity to a variation of the a width 
is constantly unity at all energies. 




Ecm (MeV) 



FIG. 7. Sensitivity of the 127 I(a,n) 130 Cs reaction cross sec- 
tions to a variation of the 7, neutron, and a width as function 
of cm. energy. The sensitivity to a variation of the a width 
is constantly unity at all energies above the threshold. 



TABLE IV. Stellar reactivity N A (cry)* for 127 I(q,7) 131 Cs as 
function of plasma temperature T. 



T 


Reactivity 


(GK) 


(cm s 


mole ) 


0.50 


1.32 x 


10 -41 


0.60 


3.79 x 


10" 37 


0.70 


1.26 x 


10" 33 


0.80 


9.34 x 


10- 31 


0.90 


2.36 x 


1Q -28 


1.00 


2.64 x 


1Q -26 


1.50 


3.04 x 


10 -19 


2.00 


5.83 x 


10- 15 


2.50 


5.00 x 


io~ 12 


3.00 


6.69 x 


10 -io 


3.50 


2.58 x 


10" 8 


4.00 


4.18 x 


10- 7 


4.50 


3.66 x 


10" 6 


5.00 


2.03 x 


10~ 5 


6.00 


2.51 x 


10- 4 


7.00 


1.41 x 


10~ 3 


8.00 


4.86 x 


10" 3 


9.00 


1.19 x 


10" 2 


10.00 


2.23 x 


10~ 2 



TABLE V. REACLIB parameters [46 
the reactivities shown in Table [TV] 



obtained from fitting 



Parameter 



(«,7) 



(7>a) 



ao 

<M 
a 2 

a 4 



-2.092540 x 10 J -2.066396 x 10 J 

-4.266646 x 10 2 -4.440888 x 10 2 

7.026704 x 10 3 
-4.701040 x 10 3 
1.388502 x 10 2 
-4.682267 



3.399136 x 10 3 



3.400636 x 10 3 



100 



discrete excited states in 130 Cs used in the H-F calcula- 
tions shown in Figs. [3] and HI Only five discrete levels, 
including the g.s., were included in the NON-SMOKER 
calculation, as given in [Q. Above the fourth excited 
state, a theoretical nuclear level density [5l[ was used. 
Version 0.8.4s of the SMARAGD code made use of an 
updated level scheme from [49) . also included in the 2010 
version of NuDAT [48j . Above the last included discrete 
state, a refitted version of the level density of [5l| is ap- 
plied and additionally an energy- dependent parity distri- 
bution [52l |. 

Using the full list of 30 levels given up to 4.04 MeV, 
the curve marked " SMARAGD (30 exc. stat.)" in Figs. [3j 
2] is obtained. Since 127 I(a,7) 131 Cs is also sensitive to 
the neutron width (see Fig. |T5]), it is also affected. The 
artificial break at about 11.8 MeV originates from the 
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FIG. 8. Ratio SMARAGD/NON-SMOKER of the newly cal- 
culated SMARAGD stellar reaction rate for 127 I(q,7) 131 Cs 
and the NON-SMOKER predictions 0. The horizontal line 
at unity ratio is drawn to guide the eye. 
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TABLE VI. Excited states in 


130 Cs used in 


the H-F calcula- 


tions. 






NON-bMOKER 144. 461 


SMARAGD [27, 45] 


Energy 


Energy 


F 


i ivie v i 


(MeV) 




a aaaa i + 
U.UUUU 1 


0.0000 


1+ 


A AO A A 0~t~ 

0.0804 2 


0.0804 


2+ 


n loi r Q + 

O.toto 2 


0.1315 


2+ 


U.14S4 I 


0.1484 


2" 


A i a o o k — 
0.1633 


0.1633 


5" 




0.2788 


6" 




0.3756 


7" 




0.5655 


8" 




0.6181 


8" 




0.8783 


9" 




0.9748 


9+ 




0.9972 


9" 




1.1265 


10+ 




1.1721 


10" 




1.2654 


10" 




1.4798 


11+ 




1.5129 


11" 




1.6738 


11+ 




1.7700 


12+ 




1.9607 


12" 




2.0748 


12+ 




2.0861 


12" 




2.1870 


13+ 




2.4468 


13+ 




2.6135 


14+ 




2.7966 


14+ 




3.0825 


15+ 




3.2498 


15+ 




3.5475 


16+ 




4.0405 


17+ 



below the (a,n) threshold and therefore there is also no 
impact on the astrophysical rate for the (a, 7) reaction. 

IV. SUMMARY AND CONCLUSIONS 

The reactions 127 I(a,7) 131 Cs and 127 I(a,n) 130 Cs have 
been studied using the activation technique combining 
X- and 7-ray countings. For the first time, measure- 
ments have been performed for these reactions close to 
the astrophysically relevant energy region, at energies 
9.50 < E c . m . < 15.15 MeV and 9.62 < E c . m . < 15.15 
MeV, respectively. Furthermore, the relative intensity of 
the 536.1 keV 7 transition was measured precisely, its 
uncertainty is reduced from 13% to 4%. 

The (a,n) data confirm the H-F predictions of the av- 
eraged a width which is essential to derive the astro- 
physical 127 I(a,7) and 131 Cs(7,a) rates, important in the 
nucleosynthesis of heavy p nuclei in a 7 process. The 
(a, 7) data reveal deficiencies in the prediction of the 7 
and/or neutron width within the measured energy range. 
This is inconsequential, however, for the stellar (0,7) and 
(7, a) rates. Despite an unchanged optical a+nucleus po- 
tential, a recalculation of the reaction rates with an im- 
proved code yielded increased rates compared to previous 
calculations at 7 process temperatures. 

There is an indication of an onset of change in the en- 
ergy dependence of the a width in the (a,n) data at the 
two lowest energies. It would be ironic if there is a reduc- 
tion of the a width at very low a energies bringing the 
new rate down to the previous value. The current data, 
however, is not sufficient to arrive at a firm conclusion. 
More data at lower energies, preferably below the (a,n) 
threshold, would be required to clarify this issue. 



inclusion of transitions calculated with a theoretical nu- 
clear level density above this energy. It is a clear sign 
that the number of included discrete levels close to this 
energy is quite different from the predicted number. This 
is due to an incomplete level scheme. As can be seen in 
Table IVII and in Ref . [49J , the given states at higher en- 
ergy are only states with higher spin from four rotational 
bands, low-spin states are not identified. Truncating the 
list of discrete levels after the fifth one yields the curve 
labeled " SMARAGD (5 exc. stat.)" in Figs. HE] It is to 
be noted that the choice of levels obviously has no impact 
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